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INVARIANCE PRINCIPLE FOR THE RANDOM CONDUCTANCE 
MODEL WITH UNBOUNDED CONDUCTANCES 

By M. T. Barlow 1 and J.-D. Deuschel 2 

University of British Columbia and Technische Universitat Berlin 

We study a continuous time random walk X in an environment of 
i.i.d. random conductances fi e £ [1, oo). We obtain heat kernel bounds 
and prove a quenched invariance principle for X. This holds even 
when E,u e — oo. 

1. Introduction. We consider the Euclidean lattice Z d with d>2. Let 
Ed, the set of nonoriented nearest neighbour bonds, and, writing e = {x, y} G 
Ed, let (n e , e G Ed) be nonnegative r.v., defined on a probability space (Q, F). 
Throughout this paper we will assume that (/x e ) is stationary and ergodic, 
and that its law is invariant under symmetries of Z d . We write \x xy = H{ x ,y} = 
ix yx , and let fj, xy = if x ^ y. Set 

(1-1) Vx = ^2^xy, P(x,y) = — . 

y ^ 

There are two natural continuous time random walks associated with \i. 
Both jump according to the transitions P(x, y). The first (the constant speed 
random walk or CSRW) waits at x for an exponential time with mean 1 
while the second (the variable speed random walk or VSRW) waits at x for 
an exponential time with mean 1/ fi x . Write Cc and Ly for their generators, 
given by 

(1-2) Ccf(x) = /i" 1 ^y(f(v) ~ /(*)), 

v 

(1.3) C v f{x) = Y,»xy{f{y)-f{x)). 
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Set 

£(f,g) = \ E E M/^) " f(v)M*) ~ a(y)). 

x&. d y&L d 

Let !/j = l,ie It is easy to check that if /, g have finite support, then 
(1.4) £(f,g) = -^9(x)^2tixy(f(y) ~ /(*)), 

x y 

and so 

(1-5) £(f,g) = -(C v f,g)u = -(C c f,g)^ 

Thus the VSRW is the Markov process associated with the Dirichlet form 
(£,£>(£)) on L 2 {y) and has stationary measure v while the CSRW is the 
Markov process associated with the Dirichlet form (£,D(£)) on L 2 (/j,) and 
has stationary measure \i. 

Let X = (X t ,t> 0, P*, x £ Z d ) be either the CSRW or the VSRW. Write 
C for its generator, 9 for its invariant measure (so either 9 = v or 9 = fx) and 
let 

(i.6) $(*,y)= PS{x : =y) 

Vy 

be the transition density of X (or heat kernel associated with C). This model, 
of a reversible (or symmetric) random walk in a random environment, is often 
called the random conductance model or RCM, particularly in the special 
case when (/x e ) are i.i.d. We are interested in the long-range behavior of 
X and, in particular, in obtaining heat kernel bounds for qf{x,y) and a 
quenched or P-a.s. invariance principle for X. When E/i e < oo, an averaged 
invariance principle is obtained in [17]. 

We begin by discussing the case when (// e ) are i.i.d. If /i e = then X 
never jumps across e. So if p + = P(/U e > 0) is less than p c (d), the critical 
probability for bond percolation in Z d , then X is P x P^-a.s. confined to 
a finite set. Thus we restrict to the case p+ > p c . A number of different 
authors have studied this model under various restrictions on the support 
of /U e . If fi e £ {0, 1} then this problem reduces to that of a random walk 
on (supercritical) percolation clusters (see [1] for heat kernel bounds, and 
[10, 29, 34] for quenched invariance principles). More generally it is useful 
to consider the following special cases: 

Case 0. c _1 < fi e < c for some c > 1; 

Case 1. < fj, e < 1; 

Case 2. 1 < n e < oo. 

For case 0, heat kernel bounds follow from the results in [18, 19], and a 
quenched invariance principle is proved in [34]. Case 1 is treated in [11, 12, 
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30]. (The papers [11, 12] consider a discrete time random walk.) These papers 
prove an invariance principle, with a strictly positive diffusion constant a 2 . 
Further, [11] shows that Gaussian heat kernel bounds do not hold in general 
in this case. 

In this paper we will look at case 2. There is not a great difference between 
the CSRW and VSRW in case 1, but in case 2, and in particular when 
E/i e = oo, the VSRW and CRSW do have different behaviors. Also, while 
the discrete time random walk with jump probabilities P(x, y) given by (1.1) 
behaves in a similar fashion to the CSRW, there is no simple discrete time 
analogue of the VSRW in case 2. We remark that our result for the CSRW 
also gives an invariance principle for the discrete time random walk with 
jump probabilities P(x,y). 

Let 

(1.7) X[ £) =eX t/£ 2, t>0. 

Our first main result is the following quenched functional central limit the- 
orem (QFCLT): 

Theorem 1.1. Let d>2. Suppose that (/j, e ) are i.i.d., and fi e > 1 P-a.s. 

(a) Let X be the VSRW. Then P-a.s. X^ £ ' converges (under P®) in law 
to a Brownian motion on M. d with covariance matrix OyL where ay > is 
nonrandom. 

(b) Let X be the CSRW. Then P-a.s. X^ converges (under P®) in law 
to a Brownian motion on M. d with covariance matrix OqL where 

2 = f a v /{2dEfi e ), ifEn e < oo, 
c \0, i/E/i e = oo. 

We also have heat kernel bounds for the VSRW. 

Theorem 1.2. Let d>2. Suppose that (/x e ) are i.i.d. and [i e > 1 P-a.s. 
Let qt(x,y) be the heat kernel for the VSRW. Let r] G (0, 1). There exist r.v. 
U x ,x S Z d , such that 

(1.8) P(£4H > n) < ci exp(-c 2 n v ) 

and constants Cj ( depending on d and the distribution of fi e ) such that the 
following hold. 

(a) For all x, y, t 

(1.9) q?(x,y)<c 3 t- d / 2 . 

(b) If\x-y\ Vt 1 / 2 > U x , then 

(1.10) q?(x,y)<c 3 t- d / 2 e- Cilx - y \ 2/t whent>\x-y\, 
qf(x, y) < c 3 exp(-c 4 |x - y|(l V \og{\x - y\/t))) 

(1.11) 

when t < \x — y\. 
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(c) Let x,yeZ d , t> 0. Then 

(1.12) q?(x,y)>c 5 t~ d / 2 e- C6lx - yl2/t if t > U 2 V \x - y\ 1+r ^ . 

(d) Let x, y£Z d and t > c 7 V \x - y\ 1+r > . Then 

(1.13) c 8 r d/2 e- C9|x - y|2/ * < Eqf{x,y) < c 10 t~ d / 2 e~ Clllx ~ yl2/t . 

Using Theorems 1.1 and 1.2 we can obtain a parabolic Harnack inequality 
(PHI) for qf (for the VSRW) by the same arguments as in [3] (see Theorem 
4.7). Since the CSRW is a time change of the VSRW, harmonic functions and 
Green's functions are the same for both processes. The PHI for the VSRW 
implies an elliptic Harnack inequality (see Corollary 4.8) which therefore 
holds for both CSRW and VSRW. Combining the invariance principle and 
the PHI, we obtain, using the methods of [3], a local limit theorem for the 
VSRW (see Theorem 5.14). 

When d > 3 the calculations in Section 6 of [3] then give bounds on the 
Green's function g u] (x,y) defined by 

POD 

(1.14) 9u{x,y)= / tf(x,y)dt. 

Jo 



Theorem 1.3. Let d > 3, and suppose that (/i e ) are i.i.d. and fi e > 1 
F-a.s. 

(a) There exist constants c±, . . . ,04 and r.v. U x , x S Z rf such that 
(1.15) F(U X > n) < ci exp(-c 2 n 1/3 ) 

and 

( 1 - 16 ) i C3 \d i - 9u>(x,y) < | C4 , ■ 2 if \x - y\ > U x A U y . 

\x — y\ \ x ~ y\ 

(b) Let C = T(f - l)/(2vr d / 2 cr^) . For any e > there exists M = M(e, u) 
with F(M < oo) = 1 such that 

(1.17) ^^<^(0,*)<^# for\x\>M(u>). 

(c) We have, F-a.s., 

(1.18) lim \x\ 2 ~ d guJ (0,x) = lim \x\ 2 ~ d E gu] (0, x) = C. 

\x\^oo \ x \~ 

Remark 1.4. (b) and (c) in Theorem 1.3 use the QFCLT, which in turn 
uses the ergodic theorem. As we do not have any rate of convergence in the 
QFCLT this means that [unlike the r.v. U x in (a)] we have no control on the 
tail of the r.v. M in (b). 
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The main difference between the RCM and the percolation case is the 
possibility of traps. Suppose e = {x, y} is a bond with /j e = K>l, and that 
all the other bonds e' adjacent to x and y have /v — 1- Then P(x,y) ~ 
1 — c/K, so X will jump between x and y O(K) times before leaving {x,y}, 
and thus the CSRW will be trapped for a time of order K in {x, y}. However, 
for the VSRW each jump takes only 0{K~ l ), so the total time spent in {x, y} 
is only 0(1). A similar effect will arise from finite clusters of bonds of high 
conductivity. 

The presence of traps of this kind is why we have, when E/i e = oo, that the 
diffusion constant a^, for the CSRW is zero. In this case it is natural to ask 
if a different scaling will give a nontrivial limit. There is a connection here 
with "aging" (see [6-8]), and in [2] it is proved that if the tail distribution 
P(/i e > t) ~ t~ a , then e a X t / £ 2 converges to the "fractional kinetic" motion 
with parameter a (see [8]). 

While we have written this paper for the case of i.i.d. conductances fj, e , 
our arguments do not require the full strength of this. In case 0, when the 
conductances are bounded (both above and below), then uniform upper and 
lower Gaussian heat kernel estimates, as in Theorem 1.2, are well known (see 
[18]). It follows (see Remark 6.3) that Theorem 1.1. holds for any stationary 
ergodic environment. On the other hand, in the unbounded case 2, there 
exist stationary ergodic environments such that the VSRW can explode in 
finite time; for an example, see Remark 6.6 below. 

For the Gaussian bounds in Theorem 1.2 we need to control the sizes of the 
clusters of high conductivity which is done by comparison between the graph 
metric d(x,y) and a new metric d (introduced by Davies in [15, 16]) which 

~ — 1/2 

is adapted to the structure of the random walk and satisfies d(x,y) < fJ, X y 
when x ~ y. This new metric is constructed by a first passage percolation 
procedure, and in this paper we have used first passage percolation argu- 
ments to compare the two metrics (see Lemma 4.2). These arguments use 
estimates from [25] which in turn use the fact that (M e are i.i.d. However, 
we could also have used a direct argument as in [12], Lemma 3.1 or [30], 
Lemma 5.3. Once we have the Gaussian bounds (with sufficiently good con- 
trol on the tails of the r.v. U x in Theorem 1.2), the quenched invariance 
principle follows with no further hypotheses on {^ e >e € E^} other than that 
it is stationary, symmetric and ergodic. Theorem 6.1 summarizes the general 
situation. 

The structure of this paper is as follows. In Sections 2 and 3 we study a 
deterministic graph T = (G,E) with edge weights fj, xy . Under certain con- 
ditions (which are P-a.s. satisfied by the VSRW on the i.i.d. RCM) we ob- 
tain heat kernel bounds in this setup. Our approach uses similar methods 
to those used in [1] for percolation clusters. However, in [1] the Carne- 
Varopoulos "long-range" bounds played an essential role at various points. 
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These bounds do not hold for the VSRW, and instead we use more general 
upper bounds obtained by Davies [15, 16], which are in terms of the metric 
d(x,y). The same metric d is also needed to control P*(d{Xt,x) > Ai 1 / 2 ) 
which is the key step in obtaining general Gaussian upper bounds. Similar 
bounds, in the context of weighted Laplacians on manifolds, are obtained by 
Grigor'yan [24]; here the metric d is the Riemannian metric. Very recently, 
and independently, Mourrat [31] has obtained upper bounds for the VSRW 
which in certain cases improve on those in Theorem 1.2. 

Once one has upper bounds, lower bounds follow by the same arguments 
as in [1], Section 5 (see Section 3). In Section 4 we then prove Theorem 1.2. 

Section 5 proves the invariance principle. We begin with the VSRW. The 
basic technique in the proof (as in many previous papers such as [10, 12, 
17, 27-29]) is to associate with Xt a process Zt on ft = [1, oo] E which is the 
environment seen from the random walk. More precisely, for each x £ Z rf , let 
T x : ft — > ft be given by 

T x (co)(z, w) = ijj(z + x, w + x). 
Assuming that Xq = we define 
(1-19) Z t {uj) = T XtH (Lo). 

One seeks to use the process Z to construct the "corrector," that is, a map 
X : ft x Z d -)• R d such that 

(1.20) M t (oj)=X t (oj)-x(uj,X t (Lu)), t>0, 

is a P^-martingale. Once one has constructed the corrector, showing the 
invariance principle for the rescaled martingale eM t j e 2 is standard, and using 
results from [12, 34], the heat kernel estimates in Theorem 1.2, together with 
the sub linear growth of %, imply that 

(1.21) lim ex(u, X t / e 2) = in P°-probability. 

However, the standard construction of the corrector is based on L 2 (P) 
calculus, which requires finiteness of the first moment of the conductance 
(see [17], page 816). In our case we wish to handle the case when E,fi e = oo, 
and so we need an alternative approach. (We remark that if E// e = oo then 
it is not easy to find suitable function spaces on ft which give a core for the 
Dirichlet form associated with Z.) 

Our solution relies on discretization. We define X n = X n , n £ Z + , and 
consider the process 

(1.22) xi e) =eX lt/£2] . 

We can control sup t<T \xj: £ ^ — x[ e \ (see Lemma 4.12), so an invariance prin- 
ciple for X^ will follow from one for X( £ ). The process X does not have 
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bounded jumps — in fact it jumps anywhere in Z d with positive probability. 
However, the long-range bounds on qf(x,y) in (1.11) give good control on 
these jumps, and, in particular, the bounds on qf(x,y) imply that 

(1.23) E££|Xi| 2 <oo, 

which is the key L 2 condition on X for the construction of the corrector. As 
we will see in Section 5, looking at the discrete time process does actually 
introduce some simplifications in the construction of the corrector x- I n the 
end (see Remark 5.15) it will turn out that the "discrete time" corrector x 
also satisfies (1.20). 

Finally, a short Section 6 makes some remarks on more general environ- 
ments, and gives an example (a one-sided spanning tree) where the process 
X fails to be conservative, and so the invariance principle fails. 

2. Transition density upper bounds on a fixed graph. Let T = (G, E) be 

an infinite (deterministic) graph, \i e) e € E, be bond conductances and v be 
a measure on G. We make the following assumptions on (G,E), /i and v. 

Assumption 2.1. (1) T is connected. 

(2) The vertex degree is uniformly bounded, 

(2.1) \{y:y~ x }\ < C D for all x £ G. 

(3) n e > for all eeE. 

(4) There exists Cm > 1 such that 

(2.2) C M l < v x = v{{x}) < C M for all x£G. 

The results of this section do not explicitly require a strictly positive 
lower bound on fj, e ; however, a later assumption [see Assumption 2.6(2)] will 
impose some control on the edges e with \i e small. 

We write [i xy for fJ-{ x , y }^ and set fx xy = if x ^ y. Let d(x,y) be the usual 
graph distance on G, and write 

(2.3) B(x,r) = {y:d(x,y)<r}. 

Let Ca < oo. We now construct, by a first passage percolation procedure, a 
second metric d on G satisfying 

(2.4) (C^ 2 V (i yy ,)\d(x,y) -d{x,y')\ 2 <1 for every x G G, y ~ y'. 

We write B(x,r) = {y:d(x,y) < r} for balls in the metric d. (In this paper 
we can take Ca = l,Jaut for possible future extensions we treat the general 
case.) To construct d we assign waiting times 

(2.5) t(e) = C A A^e 1/2 , e£E, 
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and then take d(x, y) to be the shortest journey time between x and y. More 
precisely, 

(2.6) di> )y ) = infj^( ei )J, 

where the infimum is taken over paths (ei, . . . , e n ) from x to y. Since we do 
not have a strictly positive uniform lower bound on t(e), in general there 
may not be a minimizing path. However, such paths will a.s. exist when t(ei) 
are i.i.d. positive random variables. 

Lemma 2.2. The metric d constructed above satisfies (2.4)- 

Proof. Let e = {y, z}; then d(x, z) < d(x, y) + t(e), and using (2.5) gives 
(2.4). □ 

Recall that 

W,9) = 5 £5>*y (/(l/) " f(x))(g(y)-g(x)). 

x y^x 

Let fi x = Yly^x ^xy, and extend fj, to a measure on G. Then £(f,g) is defined 
for f,geL 2 (G,fi). 

Let X = (Xt,t E [0,oo),P x ,x £ G) be the continuous time Markov chain 
on G with generator 

Cf(x)=U- 1 Y,^xy(f(y)-f(x)). 

y 

At this point we cannot exclude the possibility that X may explode, and we 
write C for the lifetime of X. Let 

II/I& =£(/,/) + 

and J- be the closure of the set of functions on G with finite support with 
respect to \\f\\e 1 ■ Then X is the Markov process associated with the Dirichlet 
form (£, F) on L 2 (G, v) (see [22]). Let qt(x, y) be the transition density (heat 
kernel) of X with respect to v. 

P x (X t = y) 

qt{x,y) = . 

Uy 

We begin by using the results of Davies [15, 16] to obtain long-range 
bounds on qt- By Proposition 5 of [15], we have 

(2.7) qt{x,y)<{y x u y T 1/2 inf exp(^(y) - ^(a) + A(^)t), 
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where A(^) = sup x b(ip,x), and 

(2.8) b{i;,x) = ^Ys^xyie^-^ +e^ v) -^ x) -2). 

x y~x 

Theorem 2.3. Assume (G,E) and /i satisfy Assumption 2.1. There 
exist constants c\,...,c± (depending only on Ca,Cd,Cm) such that the fol- 
lowing hold. 

(a) If x,y G G and D = d(x, y) < c\t, then 

(2.9) q t (x,y)<c 2 exp{-c 3 D 2 /t). 

(b) If x,y G G and D = d(x, y) > c\t, then 

(2.10) q t (x, y) < c 2 exp(-c 4J D(l V \og(D/t))). 

Proof. Fix xq, yo G G, let t > 0, and write D = d(xQ,yo). Let A > and 

set 

ipx(x) = —X(D A d(xo,x)), b(\) = supb(ip x ,x). 

X 

Let x G G, y ~ x and write = C^ 2 V /U^, 

~— 1/2 

Then as IV'a(^) — V>A(y)| < X\i xy , and cosh(x) is increasing on [0,oo), 

(2.11) J xy < 2/i X j / (cosh(A^-j / 1 / 2 ) - 1) < 2jl xy (cosh{XJl-^ 2 ) - 1). 

Using the power series for cosh we have that the right-hand side of (2.11) is 
decreasing in Jl X y, so 

J xy <Cf(e CAX + e- CAX -2). 

Hence 

b(il>x,x) < \C M C D C- 2 (e CAX + e~ CAX - 2). 

Let f(x) = e x + e~ x - 2; then b{\) < c 7 f(C A X). Thus by (2.7), and writing 
y = C A X, 

qt(x ,y ) < Cm infexp(- AD + c 7 tf{C A X)) 
(2-12) A _ 
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So if 

F( S ) = mf(-y+(2 S )" 1 (e^ + e ^-2)), 

y>0 

then 

(2.13) qt(x ,y ) < C M exp(J^-F 

and it remains to bound F. 

We have (see [15], page 70) that 

F(s) = s-\(l + 3 s ) 1 / 2 - 1) - log(s + (1 + s 2 ) 1/2 ) 

and also that F(s) < s/2(l - s 2 /W) for s > 0. Hence, if s < 3, then F(s) < 
— s/20 while if s > e, then 

F(s) < 1 - log(2s) = - log(2s/e). 

Substituting in (2.13) completes the proof. □ 

Remark 2.4. Note that if D = ct then both (2.9) and (2.10) give a 
bound of the form ce~*/ c . 

Since \i e are not bounded above, the process X may explode. The follow- 
ing condition is enough to exclude this. 

Lemma 2.5. Suppose there exists x 6 G and 9 > such that 

(2.14) y^ exp(-6d(x,y))v y < oo. 

y&G 

Then X is conservative. 

Proof. Let £ be the lifetime of X. Then as V is connected it is easy 
to see that either py(( = oo) = 1 for all y G G, or else P y ((<t)>0 for all 
y GG, t >0. 

For n > C^ 2 let /if^ = n A A^™) be the process associated with the 
conductances fi^ n \ and ql n \x,y) be the transition density of X^ with re- 
spect to v. We have qt(x,y) = lim n _ i . 00 q[ n \x, y). Note that each A^ n ) is 
conservative, and that the bounds in Theorem 2.3 hold (with the same con- 
stants Cj) for each q^ n \ With (2.2) the condition (2.14) implies that B(x,R) 
is finite for each R > 0. 

Let t > 0. With constants Cj as in Theorem 2.3, choose r large enough so 
that r > cit and c 4 (l V log(r/t)) > (9. So, if R > r, using (2.10), 

(2.15) ^2 q[ n \x,y)u y < ^ c 2 exp(-9d(x,y))u y < oo. 

y£B(x,R) c yeB{x,R) c 



D 

2C A c 7 t 
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Let e > 0; then we can take R large enough so that the right-hand side of 

(2.15) is less than e. Thus, as is conservative, for all n, 

it n \ x >y) u v> 1 - £ - 

y&B{x,R) 

So, 

P x ((>t)> V q t {x,y)v y = lim V q^\x, y)u y > 1 - e. 
y&B(x,R) y&B(x,R) 

Therefore -P X (C > t) = 1 for all t, proving that X is conservative. □ 

We now make further assumptions on the graph V and the conductances 
fj,. As we will see in Section 4, it is easy to check these for the random 
conductance model on Z d when ji e > 1. 

Assumption 2.6. (1) There exists d > 1 and GV < oo such that 

(2.16) i/(B(x,r)) < C v r d for all x G G,r > 1. 

(2) There exists a constant Cat such that the following Nash inequality 
holds. If /G L 1 (G, i/) n L 2 (G, v), then, writing ||/|| p for norms in L p {G,v), 

(2.17) ^a^^cjvii/n^ii/iir 47 "- 

Remark 2.7. Note that (2) above does place some restrictions on the 
edges e with /i e small. For example, taking x G G and f = l x , (2.17) gives 

y 

By [14] we have: 

Lemma 2.8. Suppose (2.17) holds. Then 

(2.18) q t (x,y)<ct- d/2 , x,y€G,t>0. 

To obtain better control of qt(x,y) when d(x,y) is large we need to com- 
pare the metrics d and d. Note first that d(x,y) < Ca A [i X y < when 
x ~ y, so 

(2.19) d(x,2/)<CVZ(x,y), x,yeG. 
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Definition 2.9. Let A > 1, r) g (0,1). Let x e G, r G [l,oo). We say 
(x,r) is X-good if B{x,n/X) C B{x,n) for all n > r, n G N. We say (x,Rq) 
is X-very good if for all i? > i?o 5 (y> r ) is A-good for all y G B(x,R), r > R v , 
r G N. Note that if (x,i?o) is X-very good then (x,i?i) is X-very good for all 
-Ri > -Ro- For x G G let 14 = V X (X) be the smallest integer such that (x, V^) 
is A-very good. 

Note. Unlike the definitions in [1], the event that (x,r) is A-good depends 
on the structure of T "at infinity." 

Lemma 2.10. Suppose (x,R) is X-good. 

(a) Ifd{x,y)>R, 

(2.20) A" 1 ^, y) < d(x, y) < C A d(x, y). 

(b) IfR'> (2R)V2(l + C A X)d(x,x'), then {x',R') is 2X-good. 

Proof, (a) The upper bound is given in (2.19). For the lower bound, let 
s = d(x,y) > R. Then y £ B(x,s), so y ^ B(x,s/X) and thus d(x,y) > s/X. 
(b) Let r = d(x, x'), r = d(x, x'), and s > R' . Then as s/2 > R, 

B(x',s/2X) C B{x,¥+s/2X) C B(x,X7+s/2). 

So, using (2.20), B(x', s/2X) C B(x', (1 + XC A )r + s/2) C B(x', s). □ 

Lemma 2.11. Let x G G, 9>0, r > 1. If (x,r) is X-good, then 

foou t a7( w <r J 4X)r d e- cre , ifr9>l, 

(2.21) eM-0d(x,y))u y <n^ ' - 

In particular, X is conservative. 

Proof. Write I for the left-hand side of (2.21), and D n = B(x,2 n r) - 
B{x,2 n ^ 1 r). Then 

oo oo 

(2.22) /<J^ Yl exp(-ed(x,y)/X)u y <Y,C V {^ n r) d exp(-2 n ~ 1 r9/X). 

n=ly£D n n=l 

If a > 0, d > 1 then there exists ci = c\(d) such that 



oo 



^2 2 na 

n=l 



nd„-a2" ^ J cie a , if a > 1, 
cia~ d , if a < 1, 



and using these bounds in (2.22) completes the proof. □ 
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Lemma 2.12. Let x G G and suppose that (x,r) is X-good. If t € (0,1), 
then 

(2.23) E x d(x,X t ) p <c(X)r p+d . 

Proof. Using the bound (2.10) a similar calculation to that in Lemma 
2.11 gives 

oo 

E x d(x,X t ) p < r p + J2C V {2 n r) d+p e- cd{x ' y) 
n=l 

oo 

< r p + cr d+p 2 ™ (d+p) exp(-c'2 n r/A) < cr d+p . 

n=l 

We now follow the arguments in [1], Section 3 (the "Bass-Nash method") 
to obtain Gaussian upper bounds on qt(x, y). As in Lemma 2.5 for 1 < n < oo, 

let \ily = fJ-xy A n, be the associated VSRW, and q( nS, (x,y) be the 

transition density of X^ n \ Let xq E 67, and set 

(2.24) M n (t) = M n (x Q ,t) = E Xo d(x ,xi n } ) =Y,d(x ,y)q < f l) (x ,y)v y , 

y 

(2.25) Q n {t) =Q n (x ,t) = - ^2q^ n) (x ,y)logql n \x ,y)uy. 

y 

The following three inequalities lead, by Nash's argument [32], to upper 
bounds on M n (t) which are uniform in n. [We remark that we only need 
the approximations X^ to justify an interchange of sums in part (c) of the 
following lemma.] 

Lemma 2.13. Let xq^G and r > 1. Suppose (xo,r) is X-good, and 1 < 
n < oo. There exist constants Ci, independent of n, such that the following 
hold. 

(a) We have, for t>0, 

(2.26) Q n {xQ,t)> Cl + \d\ogt. 
(b) 

(2.27) M n (x Q ,t) > C2e Q^o,t)/d provided either M n (x ,t) >r ort> c 3 r 2 . 
(c) For t > 0, 

(2.28) Q' n (t)>c A M' n (t) 2 . 
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Proof. We write Q n (t), M n (t) for Q n (xo,t), M n (x ,t). (a) is immediate 
from (2.18) and the fact that since is conservative, ^2 y qt n \xo,y)^ y = 1. 

(b) The proof is similar to those in [1, 4, 32]. First note that (2.16) and 
Lemma 2.8 give that 

(2.29) M n (t) > r provided t > c 3 r 2 . 

By Lemma 2.11, provided ar < 1, 

y£G y£B(x ,r) y<£B(x ,r) 

< cr d + ca~ d < ca~ d . 

Now uilogu + 9) > —e~ x ~ e for u > 0. So, setting = ad(xQ,y) + 6, where 
a < 1/r, 

-Qn{t) +aM n (t) + b= ^2 q[ n) (x ,y)(logq ( t n \x ,y) + ad(x ,y) + b)u y 

> _ e -l-6^ g-ad^o*)^ > _ C5e -V d . 
yGG 

Setting a = 1/M n (i) and e 6 = M n (i) d = o" d , we obtain 

-Qn(0 + 1 + d\ogM n {t) > -ce 
and rearranging gives (b). 

(c) Set = g t (n) (x ,a;), and let b t (x,y) = f t (x) + / t (y). We have, using 
(2-4), 

M' n {t) = ^2d(x ,y) ^^-vy =J2d(x ,y)C^f t (y)u y 
v y 



~^2^2fJ>xy(d{x ,y) - d(x ,x))(f t (y) - f t (x)); 



2 

x y 



the final interchange of sums can be justified using (2.9) and (2.10) and the 
fact that is uniformly bounded. So using (2.4), 

Kit) <^E(^l d >o,1/) - d(«o,x)|^(x,y) 1 / a )(^ l/ «ff "^ a)l 



x yr^x 



a EE(M ^)-,(^M|i_^ 

x y^x 
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Y /2 ^^ m) - ftW) 2 \ 1/2 



^ \ (E Z bt ( x ' f )) (Z Z ^ - 



Now 

(2.30) X) Z) ^ f ) = 2 Z) Z) ^ ^ 2C D CW Z M*)"* = 2C dCm- 

x y^x x yr^x x 

So, 

M'ft)kcVVii (/ f (|/)-/ f (x)) 2 

Since we have, for u, v > 0, 

< (u — v) (log u — log v) , 

we deduce 

M^) 2 < J>^(/ t (y) - f t (x))Qogf t (y) - log/ t (x)). 

X J/ 

Thus 

Qn(*) = - Z^ 1 + lo s Mv))^ n) Mv>v 

(2.31) = i X)E^( lo s/*(f) -hsMxMMv) -ft(x)) 

x y 

>cM^t) 2 ; 

where again the interchange of sums uses (2.9) and (2.10) and the fact that 
/i^ is uniformly bounded. □ 

Proposition 2.14. Let x e G, r > 1 and (xo,r) be X-good. Then 

(2.32) cit 1/2 <M oo (x ,*) <c 2 t 1/2 fort>c 3 r 2 . 

Proof. Note that the lower bound is immediate from Lemma 2.8. For 
the upper bound let first n < oo. Set R n {t) = d~ 1 (Q n (t) — c\ — ^dlogt), so 
that R n (t) > by (2.26). Then 

M n (t) = J M' n {s)da<cJ Q' n (s) 1/2 ds<c J (K( s ) + ^j ' ds - 
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Using the inequality (a + fe) 1 / 2 < ft 1 / 2 + a/ {2b) 1 / 2 and integrating by parts 
we obtain 

M n (t) < ct 1 / 2 + c [ s l,2 R' n {s)ds 
Jo 

< ct 1 ' 2 + c(l + R n (t)t 1 / 2 ) < ct l / 2 (l + R n (t)). 

By (2.27) we also have M n (t) > t l / 2 e Rn( ^ for t > cr 2 . Thus R n (t) is bounded 
for t > cr 2 and this implies that 

(2.33) dt 1/2 < M n (x Q ,t) < c 2 t 1/2 for t>c 3 r 2 . 

Since the constants in Lemma 2.13 are independent of n, the constants q 
in (2.33) are also independent of n. Since M^it) < lim inf n ^oo M n (t), (2.32) 
then follows. □ 

Lemma 2.15. Let x £ G, r > 1 and (x,r) be X-good. Then 

(2.34) cit l/2 < E x d(x, X t ) < c 2 t 1/2 for t > c 3 r 2 . 

Proof. Since d(x,X t ) > C^ 1 d(x,X t ), the first inequality is clear. Let 
C3 be as in Proposition 2.14, and t = C3R 2 , so R\ > r. Then if A = B(x, Ri) c , 
using (2.20), 

E x d(x,X t ) < \R l + E x {d(x,X t )-X t e A) 
<\R 1 + E x (Xd(x,X t );X t £A) 
< XRt + Xc 4 t 1/2 < crj 1 / 2 . □ 

The next few results follow quite closely along the lines of [1]. Let 
r(x, r) = inf{t : d(x, X t ) > r}. 



Lemma 2.16. There exist constants c\, c 2 , C3 such that if R>c\ and 

(2.35) (y,C2-R) is X-good for all y S B(x,R), 
then ift = R 2 /(2c 3 ) 

(2.36) P x (r(x,R/2)<t )<± 
and hence for t>0, 

(2.37) p*(rfa R )<t)<l + rt. 
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Proof. Write r = t(x,R/2), and c[ for the constants q in Lemma 2.15. 
Let C3 = 64(c' 2 ) 2 . Choose c 2 so that r = c 2 R satisfies c' 3 r 2 = to, and c\ so that 
r > 1. Then as (2.35) holds we can use Lemma 2.15 to bound E y d{y,X s ) for 
s>to,y£ B(x, R). So, 

c' 2 (2t ) 1 / 2 > i^(x,X 2io ) > X toAr ) -d(X toAr ,X 2to )) 

>£; x 'l (T < io) (i(x,X T ) - sup sup E y d(y,X 2to - s ) 

y£B(x,R) 0<s<t 

>P x (T<t )R/2-c' 2 (2t Q ) l l 2 , 

and rearranging we obtain (2.36). 

Inequality (2.37) now follows easily; if t < t Q , P x (t(x, R)<t)< P x (t(x, R/ 
2) < to) < ^ while if t > to, then the right-hand side of (2.37) is greater than 
1. □ 

To obtain the Gaussian upper bound on qt(x,y) we need to prove that 
the process X does not move too rapidly across a ball B(x,R). We choose 
r <C R, and use the fact that if X moves across B(x,R) in the time interval 
[0, t], then it has to move across many smaller balls of side r in the same 
period; the estimate (2.37) is enough to bound the probability of this. Our 
argument uses the following easily proved estimate: 

Lemma 2.17 (See [5], Lemma 1.1). Let £i,£2> • • • >£n> ^ ^ e nonnegative 
r.v. such that V > Yli£i- Suppose that for some pE (0,1), a > 0, 

(2.38) P(&<t|<7(£i J ... ) &_i))<p + <rf J t>0. 
Then 

( ant \ 1/2 

(2.39) \ogP{V <t)<2[^\ -nlog(l/p). 

Proposition 2.18. There exists constants c\, .. .,04 such that if x G G, 
R>c\, t> c\R and 

(2.40) (z,c 2 t/R) is A- good for all z£B(x,R), 
then 

(2.41) P x {t(x, R)<t)< c 3 e~ C4R2/t . 

Proof. Let 1 < m < R/2, and set r = R/2m. Define stopping times 

S = 0, S i =M{t>S i - 1 :d(X Si _ 1 ,X t )>r}, i>l. 

Set Tj = Si — Si-i, and write Ft = cr(X s ,s < t) for the filtration of X. As 
d(X Si ,X Si+1 ) < r + 1 < 2r, we have S m < t(x,R) and X Si G B(x,R) for 
< i < m — 1. 
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Suppose for the moment that m is such that we can apply Lemma 2.16 
to control each n. Then 

(2.42) p»( T .< u |jr )< + 5 u >0,l<i<m, 

2 r z 

so writing p = ^, a = c^/r 2 and using (2.39), we obtain 
logics, r) < t) < log f £ r, < tj 

< 2(amt/p) 1 ^ 2 — mlogp -1 

(2.43) 

£ - C6m ( 2 -( £ f 5 ) 1/2 ) 

= -c 6 m(2- (m/K) 1/2 ), 

where re = R 2 /(cjt). If re is such that we can choose m S N with K<m< 2k, 
and so that (2.42) holds, then (2.43) implies (2.41). 

We can choose c\ so that re < R/2 — 1. If re < 1 then, adjusting the constant 
C3 appropriately, (2.41) is immediate. If 1 < re < R/2 — 1 then let m = |_rej + 
1 < 2k. Then ^c§{t/R) <r< ^c^(t/R), and so choosing c 2 suitably, (2.40) 
and Lemma 2.16 imply (2.42). □ 

Recall that V x is the smallest R such that (x,R) is A-very good. 

Theorem 2.19. Let x,y € G, and write D = d(x, y). Suppose that either 
D>C!VV X ort>D 2 . 

(a) If c 2 D <t, then 

(2.44) q t (x, y) < c 3 r d / 2 e" *^)' '/* . 

(b) // c 2 D > t, then 

(2.45) q t (x,y)<c 3 exp(-c 4 D(lV\og(D/t))). 

Proof. We need to consider various cases. First, if t> D 2 , then (2.44) 
follows from (2.18). So we can suppose D > V x . Note that by (2.20) d(x,y) > 

If t < c 2 D, then (2.45) follows from (2.10). If c 2 D<t< c 6 D 2 /logD, then 
Theorem 2.3 gives 

Qt(x,y) <c 8 exp(-2c 7 D 2 /t). 
Choosing c& small enough we have exp(-c 7 L> 2 /0 < t~ d l 2 and (2.44) follows. 



INVARIANCE PRINCIPLE FOR THE RANDOM CONDUCTANCE MODEL 19 



It remains to consider the case cqD 2 / log D <t < D 2 . Let A x = {z: d(x, z) < 
d(y,z)}, A y = G-A x , t' = t/2, D' = D/2. Note that B(x,D')QA x . Then 

(2.46) v x P x (X t = y) = v x P x (X t = y,X t , G A y ) + v x P x (X t = y,X t , G A x ). 
To bound the first term in (2.46), and writing r = t(x,D'), we have 

P x (X t = y,X tl eA y ) 

= P x (T<t',X t/ eA y ,X t = y) 

(2.47) <E x (l {T<tl} P x ^X t . T = y)) 
<P x (t(x,D') <t') sup q t -u(z,y)v y 

z£dB(x,D>),u<t/2 

<ct- d / 2 P x {T{x,D')<t/2). 
Similarly, using symmetry, for the second term in (2.46) we have 
v x P x (X t = y, X t , G A x ) = u y py(X t = x, X t , G A x ) 

(2.48) 

<ct^ 2 P y {T(y,D')<t/2). 

It remains to verify that we can use Proposition 2.18 to bound the terms 
P z (t(z, D') < t/2) for z = x,y. Writing for the constants Cj in Proposition 
2.18, taking c\ large enough we have D' > c[, and t' > d x D' . As (x, V x ) is 
A-very good, and D > V x , we have that (z,r) is A-good for z G B(x,2D), 
and r > (2D) r > . We have c' 2 t'/D' > (2D) r > provided t > c$D 1+r i , and since 
t > cqD 2 / log D this holds by adjusting the constant c\. So 

(2.49) P z (t(w,D') <t/2) <cexp(-c'D 2 /t) iorz = x,y, 
and combining the estimates (2.46)-(2.49) completes the proof. □ 

Remark 2.20. This theorem does not give any bound for qt(x,y) when 
D < c\ V V x and t < D 2 . In this case we still have the global upper bound 
(2.18). In addition the "long-range" bounds in Theorem 2.3, bound qt(x,y) 
in terms of d(x,y), but we do not have a bound in terms of D. 

The final result of this section is that, under fairly mild additional condi- 
tions, functions which are harmonic for the discrete time process X n , n G Z + 
are also harmonic for the continuous time process Xt , t G 1R+ . At the end 
of Section 5 we will use this remark to note that the corrector constructed 
using the discrete time process also gives us a corrector for the continuous 
time process X t . 

Let X be the discrete time process given by X n = X n , n G Z+. Write 

(2.50) Cf(x) =^2qx(x,y)v v (f(y) - f(x)). 

y 
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We say h is C harmonic if the sum in (2.50) converges absolutely for all x, 
and Chix) = for all x. This implies that (h(X n ),n £ Z+) is a P x -martingale 
for each x £ G. 

For x £ Cr let k x = fj, x /v x be the jump rate out of x by X. Set 

(2.51) A(i^) = {y£G:K ?/ <i^}. 

Lemma 2.21. Let V satisfy Assumption 2.6. In addition suppose that 
there exist {xq,Tq) such that (xoi r o) *s Xgood, and that there exist Rq, K 
such when R> Rq then every self avoiding path 7 from xo to B(xq,R) c 
contains at least R lj/2 points in A(K). Let h:G^-R be £ harmonic, and 
satisfy the growth condition 

(2.52) \h(x)\ <Ci + Cid(x ,x) p 

for some p £ [l,oo). Then Ch = 0, so that h is harmonic for X . 

PROOF. By Lemma 2.10(b) we have that (x,r) is 2A-good if r/2 = ro V 
(1 + CA^)d(xo,x). So by Lemma 2.12 there exists C2 (depending on C\ and 
A) so that if s £ [0, 1], 

(2.53) E x \h(X 8 )\ <cr d +P< C 2 (r p +d + d(x , xf); 

it follows that E x h(Xt) is well defined for any t>0. Set for s £ [0,oo) 

h s (x) = E x h(X s ). 

To prove the lemma, it is sufficient to prove that h = h s for every s; this 
implies that h{X t ) is a continuous time martingale and hence that Ch = 0. 
We have 

h s+1 (x) = E x {E x 'h{X l )) = E x (h{X s )) = h s (x), 

so s — > h s has period 1. We extend h s by periodicity to s £ M. Since E x h s {X\) = 
E x h{Xi +s ) = h s (x), each h s is >C-harmonic. Let 

k(x)= sup \h S2 (x) - h Sl (x)\; 

0<si<s 2 <l 

note that by (2.53) we have 

(2.54) k(x) < 2 sup E x \h(X s )\ < 2C 2 (r p+d + d(x , x) p ) . 

8<1 

Fix x £ G, and write k = k x . Write P xy = fj, xy / ' ji x for the jump probabil- 
ities of X. Then by conditioning on the time of the first jump of X, if it 
occurs in [0, 1], we obtain 

h s (x) = e~ K h s (x) +^2,P xy f Ke~ KU h s _ u (y) du. 
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So 

h 8 (x)(l 



(2.55) =IZ jP ^( K ^ e K h s -u(y)du + K(e KU - e K )h s _ u {y) duj 

= J2 P *y( K J e~ K h u (y)du + K(e~ KU - e- K )h s _ u (y) di?J . 

Then (2.55) implies that 

(2.56) k(x)(l - e-"*) <y^P xy k(y)(l - (1 + K x )e~ K *). 



y 



So if k{x) > 0, then there exists y ~ x with k{y) > k(x). Further, if k x < K, 
then there exists 5 > (depending only on Cm and ET) such that 

(2.57) k(y) > (1 + S)k(x) for some 

Suppose now that there exists x\ with k{x\) > 0. Then there exists a non- 
intersecting infinite path 71 starting at x\ on which k is strictly increasing. 
Let 72 be a shortest path from xo to a closest point y on 7 to Xq, and let 
D be the length of 72. Combining 72 and the infinite segment of 71 starting 
at y, we obtain a path 7 = (a^^i, . . .) for which k(z n ) > for all n > D. 
Let R> Rq, and let % be the first point in 7 n B(xo,R) c n ^4(-fQ- Then 
i?i = d{xQ,WR) > R. So, using (2.54), (2.57) and the condition on A(K), 

2C 2 (r p +d + i?f) > k(w R ) > (1 + «5)^ /2 - D fc( a;i ), 
which is a contradiction if i? is large enough. □ 

3. Lower bounds and Harnack inequalities. Unlike the papers [10, 12, 34] 
we will need to make explicit use of heat kernel lower bounds in our proof 
of the invariance principle Theorem 1.1 (see Lemma 5.9). 

In this section we specialize to the case when L is the d-dimensional 
Euclidean lattice, and fj, e are bond conductances with fi e >l. We continue 
to assume that Assumptions 2.1 and 2.6 hold. Note that balls and distance 
are with respect to the graph distance on Z d . 

We can follow the arguments in Section 5 of [1] fairly closely. First, as 
fj, X y > 1 when x ~ y, by comparison with the standard Dirichlet form £q on 
Z d we have a weighted Poincare inequality as in [1], Theorem 4.8. 

Theorem 3.1. Let B = B(x ,R), p B {y) = d(y,B c ) and cp(x) = R 2 (R A 
PB (y)) 2 . Then iff-.B^M, 

(3.1) inf^(/(z)-a)V(*K<C7? 2 £ {f{x)-f{y))\{x)Mp{y),i xy . 

xdB x,y£B 
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Using this, and the method of Fabes and Stroock [21] we obtain a lower 
bound of the form qt(x,y) > ct~ d / 2 when x,y are close enough together. 

Proposition 3.2. Let xq G Ij d and R>c\. Then provided 

(3.2) {z, c 2 R) is X-good for all z G B(xo,R), 
we have 

(3.3) q t (xi,x 2 )>c 1 t- d/2 for Xi,x 2 e B{x Q ,R/2),\R 2 <t< R 2 . 

Proof. This can be proved using the argument in [1], Proposition 5.1, 
with only minor changes. Note that we need to show that ¥ Xl (Xt ^ B(xq, 
2R/3)) < \ when t = OR 2 and 6 is sufficiently small (see (5.2) and (5.9) in 
[1]). [There is a missing minus sign in exponential in the last line of (5.2).] 
This is done using Lemma 2.16, and so to satisfy (2.40) we need (3.2). □ 

Theorem 3.3. Let x,y G t > 0, and write D = d(x,y). Suppose that 

(3.4) t > ci V V 2 V D 1+r >. 
Then 

(3.5) qt(x,y)>c 2 t^ 2 e- csd ^ 2 / t . 

Proof. The proof as in [1], Lemma 5.2, Theorem 5.3, follows by a stan- 
dard chaining argument. We just give the details of the conditions on V x , D 
and t needed to make this argument work. 

First, if D 2 < t then the lower bound in (3.5) is just t~ d / 2 , so we can use 
Proposition 3.2. We set R = ct 1/2 . Then t > V 2 implies R > V x , so B(x,R) 
is A-very good, and so as cR > R v , (3.2) holds. 

If D 2 > t then we set R = 2D, r = ct/D. We apply Proposition 3.2 in a 
chain of balls B{ = B(zi,r) linking x and y. (See [1], Lemma 5.2, or [21] 
for details of the calculations.) Since D > t 1 / 2 > V x , we have that (x, R) is 
A-very good, and hence that (z,cr') is A-good for all r' > R v , z G B(x,R). 
As r = ct/D > dRP, (3.2) holds for all the balls B { . □ 

Remark 3.4. 1. Note that the lower bounds do not extend to the range 
when tc±D. The difficulty is that if t ~ D, then we need Proposition 3.2 for 
a chain of balls of radius 0(1) connecting x and y. The hypothesis "very 
good" is not enough to ensure this. 

However, the chaining argument does not need (3.2) for all points in 
B(xo,R), but just for a suitable chain connecting x and y. In [1] this fact 
was used to obtain full Gaussian lower bounds. It is likely that the same 
approach will work for the random conductance model, but we do not pursue 
this point, since the bounds in Theorem 3.3 are enough for most applications. 
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2. A well-known theorem (see [23, 33]) states that for Brownian motion 
on a manifold Gaussian bounds are equivalent to two conditions: volume 
doubling plus a family of Poincare inequalities. This theorem was extended 
to graphs in [18]. Since we have volume doubling (for u) and the Poincare 
inequalities hold (since /i e > 1), one might therefore ask if Theorems 2.19 
and 3.3 follow immediately from known results. 

However, it is clear that some conditions on fi e are needed before The- 
orem 2.19 holds — one has to prevent X from moving a long distance in a 
very short time. In fact, examination of the theorems in [18, 23, 33] shows 
that in each case there is a "hidden" additional assumption which prevents 
the process from moving too quickly. For example, [18] considers a discrete 
time nearest neighbour random walk. 

For BcZ d let qj^(x, y) be the transition density for the processes X killed 
on exiting from B. 

Lemma 3.5. Let (xq,R) be X-very good. Then 

(3.6) qf [x °' R \x,y)>c 1 t~ d/2 , x, y E B(x , 3R/4), c 2 R 2 < t < R 2 . 

Proof. Using Theorem 2.19 and Proposition 3.2, this follows, as in [1], 
Lemma 5.8, by the argument in [21], Lemma 5.1. □ 

We now give a parabolic Harnack inequality (PHI) for X . The statement 
requires a little extra notation. If A C Z d we write dA = {y:y ~ x for some 
x E A} for the exterior boundary of A, and A = A U dA. We call a function 
u(t,x) caloric in a space-time region Q = A x (0,T) C [0,oo) x Z rf if u is 
defined onQ = Ix[0, T] and 
d 

— u(t, x) = Cvu(t, x), (t, x) G Q. 

Write Q(x,R,T) = B(x,R) x (0,T], Q_(x,R,T) = B(x,\R) x [±T, ±T] and 
Q+(x,R,T) = B(x,±R) x [lT,T]. 

Definition 3.6. We say the parabolic Harnack inequality (PHI) holds 
with constant Cp for Q = Q(x,R,T) if whenever u = u(t,x) is nonnegative 
and caloric on Q, then 

(3.7) sup u(t,x)<Cp inf u(t,x). 
(t,x)eQ-(x,R,T) (t,x)eQ+(x,R,T) 

Theorem 3.7 (Parabolic Harnack inequality). There exists a constant 
Cp such that if (x,R) is X-very good. Then the PHI holds with constant Cp 
in Q{x,R,R 2 ). 

Proof. Using the heat kernel bounds in Theorems 2.19 and 3.5, and 
Lemma 3.5, this follows by the same argument as in [3], Theorem 3.1. □ 
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4. Heat kernel bounds for the RCM. In this section we prove Theorem 
1.2. Let Ed be the edges of the Euclidean lattice Z d , and let Q, = [l,oo] Ed . 
Let P be a probability measure on which makes the coordinates i.i.d. 
with a law on [l,oo). We set fi e (io) = uj(e) for e G Ed, and for each ui G VL 
we consider the random walk X on the graph i^L d ,Ed) with conductances 

Using the notation of Section 2 we take u x = l for all x, so that we can 
take Cm = 1- We write P% for the law of X started at x, and 

tf(x,y) = PZ(X t = y) 

for the transition density of X . 

Lemma 4.1. The graph (Z^E^), conductances [i e and random walk X 
satisfy Assumptions 2.1 and 2.6 with ¥ '-probability 1. 

Proof. Assumption 2.1 is immediate; note we can take Cp = 2d. Set- 
ting Cy = 2 d Assumption 2.6(1) is also immediate. 

Since we have /x(e) > 1 for all edges in if £q is the Dirichlet form of 
the standard continuous time random walk on Z d , then £(f,f) > £o{f,f), 
so that the standard Nash inequality on Z rf (see [14]) implies Assumption 
2.6(2) with a constant Cm depending only on d. □ 

In what follows we set Ca = 1; thus none of the constants Cd,Ca, Cm, Cat, Cy 
depend on the law of fi e (apart from the fact that P(^i e G [l,oo)) = 1). 

Let d(x,y) be the graph metric on (Z d ,^), and d = d(u) be the metric 
given by (2.6); as in the previous sections we write B(x,r) for balls in the d 
metric. Write Be(x, r) = {y G ¥l d : \ x — y\ < r} for the Euclidean ball center 
x and radius r. 

Lemma 4.2. There exists a constant Xq > such that 

(4.1) F(B(0, r) C B E (0, A r)) > 1 - Cl e~ C2r '. 

Proof. We use results on first passage percolation from [25]. As in [25] 
let 6o,n be the first time B(0,t) reaches the hyperplane {xi = n}. Using [25], 
Theorem 2.18, there exists /io such that lim n ?i~ 1 6o,n = Mo> a - s - an d in L • 
By [25], Theorem 1.15, we have pio > 0. By [25], Theorem 5.2, there exist 
C3 , C4 > such that 

(4.2) P(6 ,n < |np ) < c 3 e~ Cin , n > 0. 

The times for B(0, t) to hit each hyperplane {xi = ±n}, for i = 1, . . . ,d have 
the Scinie law cis 6o n •, so we deduce 

(4.3) P(B(0, Ifion) C [-n, n] d ) > 1 - 2dc 3 e" C4n , n > 0, 
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and (4.1) follows easily. □ 



Note that Ao does depend on the law of fj, e . We fix 77 € (0, 1), and define 
good and very good as in Section 2, with A replaced by Ao; and we write V x 
for the smallest integer such that (x, V x ) is very good. 

Theorem 4.3. (a) 

P((x,r) is not good) < ce~ cr , r > r$. 

(b) 

(4.4) P(14>n)<cexp(-cn r '). 

Proof. Let G(y, r) = {{y, r) is good}, and F(R) = {(y, r) is good for all 
y<EB(Q,2R),r>RP}. Then 

oo 

W(G(y,r) c )<Y,ce- cn <ce- cr . 

n=r 

So, 

oo 

¥(F(R) C ) < cR d c 3 e- C2fe < cexp(-c J R"), 

k=R r i 

and since {V > n} = \J™ F(k) c , (b) follows. □ 
Using Lemma 2.11 we obtain the following: 
Corollary 4.4. X is conservative with P '-probability 1. 

Corollary 4.5. Let xeZ d . Then 

(4.5) lim lim P°(\X t \ >My/i) = 0, F-a.s. 

M— >oo t->oo 

PROOF. By Lemma 2.15, for t > cVq 2 (o;), 

P°(\X t \ >My/i) KcM-H-WEft^Xt) <cM~\ 
and (4.5) follows. □ 

Theorem 4.6. There exist r.v. U x ,x G Z d , such that 

(4.6) P(£4H > n) < a exp(-c 2 n' ? ), 



26 M. T. BARLOW AND J.-D. DEUSCHEL 

and if \x — y\ V t 1 ^ 2 > U x , then 

(4.7) <g(x,y) < c 3 t- d/2 e- C4|a; - y|2/t whent>\x-y\, 

(4.8) qf(x, y) < C3 exp(— c^\x — y\(l V log(|x — y\/t))) when t < \x — y\. 
Further, 

(4.9) q y{x 1 y)> C Qt- d l 2 e- c ^ x - y \ 2 l t if t > U 2 V \x - y\ 1+v . 

Proof. We take U x = Cs(V x + 1) where c$ > 1. The bounds then follow 
from Theorems 2.19 and 3.3. [Note that the bounds (4.7) and (4.8) are of 
the same form if d(x,y) <t< cd(x,y).] We use the constant c$ to adjust 
between the the Euclidean metric \x — y\ and the graph metric d(x,y), and 
to absorb the conditions d(x,y) > c and t > c into (4.6). □ 

Theorem 4.7. There exists a constant Cp and r.v. U x ,x G 1^ with 

(4.10) P(C/xH > n) < ci exp(-c 2 n' 7 ), 

such that if R>U X then a PHI with constant Cp holds for Q(x,R,R 2 ). 

Proof. This is immediate from Theorem 3.7 and (4.4). □ 

The PHI implies an elliptic Harnack inequality (EHI), which holds for the 
CSRW as well as the VSRW. A function h is harmonic on A C 7L d if it is 
defined on A and Lyh{x) = [or equivalently Cch{x) = 0] for x G A. 

Corollary 4.8. There exists a constant Ch and r.v. U x ,x £ Z d with 

(4.11) F(U x (oj)>n)<c 1 exp(-c 2 n 11 ), 

such that if R> U x , then an EHI with constant Ch holds for B(x,R); if 
h > is harmonic in B(x,R), then 

(4.12) sup h(y)<C H inf h(y). 

yeB(x,R/2) y&B{xfi/i) 
We have the following averaged bounds: 

Theorem 4.9. (a) Let x,y <E Z d and t > a V \x - y| 1+r? . Then 

(4.13) c 2 t- d/2 e- C3lx - yl2/t < Eq?(x, y) < c^-^V^-^/' : . 
(b) We have 

(4.14) EE°jX t \ 2 <c 6 t, t>l. 
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Proof, (a) Let D = \x — y\. Choose c\ so that F(U X > c\ /2 ) < \. Then 
if t>ci VD 1+r >, by (4.9), 

Eq?(x,y) > E( q ?(x,y);U 2 x < Cl ) > \ct~ d l 2 e'^ V* . 

For the upper bound, let 7/ = 1 — 77, and -R^. be the r.v. given in Theorem 
4.6 using 7/ instead of 77. Then by (4.7) and (4.6), 

Eqf(x,y)=E(qf(x,yy,R' x > D) + E(qf(x, y); R' x < D) 
< ct- d / 2 e- cDV ' + ct-^e-^ 2 / 1 . 

Since the second term is larger, the upper bound in (4.13) follows, 
(b) We have 

(4.15) E°jX t \ 2 = ^2\x\ 2 qf(0,x). 

X 

We split the sum in (4.15) into three parts. First, 

(4.16) Y MV(0,x) <£/<?• 

|x|<Uo 

Next, using (4.7), 

(4.17) y mVm< J2 \A 2 ti~ m z~ cWlt <ct. 

Uo<\x\<ct Uo<\x\<ct 

Finally, using (4.8), 

(4.18) Yl \x\ 2 qt{0,x)< Y k| 2 ce~ c|x| < c'. 

ctVUo<\x\ ct<\x\ 

Combining (4.15), (4.16) and (4.17) gives 

E°\X t \ 2 <ct + c'Ul 
and as by (4.6) EU$ < 00 and t > 1, we obtain (4.14). □ 

Remark 4.10. Combining Lemma 2.8, Theorems 4.6 and 4.9 completes 
the proof of Theorem 1.2. 

Now let 

(4.19) Z t (e) =eX t/e 2, 0<e<l. 
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Theorem 4.11. Let T > 0, 5 > 0, r > 0. Then 

(4.20) lim supP°(sup|X s (e) | > r) -> 0, 

(4.21) lim lim sup P°( sup \X$ - X$\ > r) = 0. 

<*-»0 >i-s 2 |<<5,Si<T 7 

PROOF. By Theorem 4.6, if R > U , then 

P°(supX s > p) < cexp(-cP 2 / T )- 

\s<T f 

So if R > U , then R/e > U and 

P°fsup|xi e) | > R) = P°( sup \X S \ > R/e) < cexp(—cR 2 /T), 

^s<T ' \<T/e 2 ' 

proving (4.20). 

To prove (4.21) write 

(4.22) p(T,5,r)=W sup \X S2 -X Sl \>r 

K \si-s 2 \<S, Si <T 

so that 

P°( sup |X| 2 £ ) > r ) = p(T/e 2 ,S/s 2 ,r/s). 

v |si-s 2 |<<5,Si<T 7 

We begin by bounding p(T,5,r) for fixed T, 5 and r. Let k £ (0, |), = 
m a x a:e-B(o,fl) ^4, and ff(-R) = {C/^ < R K }. Then 

(4.23) F(H(R) C ) < cR d exp(-cR Kri ), 

so by Borel-Cantelli there exists Rq = Ro(uj) such that oj £ H(R) for all 

R ^ Po- 
Let 

(4.24) Z k = swp \X kS+s -X k5 \. 

0<s<5 

Then if K = \ T /5\ and Z* = maxo<fc<_ft- Z k , it is enough to control Z* since 

sup \X S2 -X S1 \ < 2Z*. 

|si-S2|<5,Si<T 

Let R> 1. Then 

(4.25) P°(Z* > r) < P°(r(0,P) < T) + P°(Z* > r,r(0,P) > T). 
By Proposition 2.18 we have 

(4.26) P°{t(0,R) <T)< ceM~cR 2 /T), 
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provided that (0, R) is very good. For this it is sufficient that R > Uq(lo). 
Now, 

K 

P°(Z* > r,r(0,R) > T) < J]P°(Z fc > r,X kS € B(0,R)) 

k=0 
K 

<E E P%(T(y,r)<5)P°(X kS = y). 

k=0 y eB(0,R) 

Again by Proposition 2.18, for y £ B(0,R), 

(4.27) P£{t{v,t) <S)< ceM~cr 2 /S), 

provided r>U* R . This will hold if R > Rq(lo) and r> R K . Combining (4.25), 
(4.26), (4.27), we obtain 

(4.28) P°(Z* >r)< cexp(-cR 2 /T) + c(T/S) exp(-cr 2 /5), 

provided R>Ro(uj) and r>R K . 
Hence 

(4.29) p(T/e 2 ,5/e 2 ,2r/e) < cexp(-ci? 2 /T) + c(T/S)exp{-cr 2 /5), 

provided R > eRq and r > R K e l ~ K . For fixed r, 5 choose R so that R> Rq 
and R 2 /T > r 2 /5. Then 

p{T/e 2 ,5/e, 2r/e) < cT5~ l exp(-cr 2 /<5) when e l ~ K < rR K . 

Hence 

lim supp(T/e 2 ,S/e,2r /e) < cTcT 1 exp(-cr 2 /<5) , 
and (4.21) follows. □ 

For n G N let X n = X n , and set 

(4.30) xi £) =sX We 2 } , 0<e<l. 



Lemma 4.12. Foranyu>0, 
(4.31) lim P° ( sup \X& - \ > u) = 0. 

e^O V < s < r / 

Proof. In the notation of the previous theorem, it is sufficient to bound 
p(T/e 2 ,l,u/e); using (4.29) we have 

p{T/e 2 , 1, u/e) < cTe~ 2 exp(-cu 2 /e 2 ), 

provided there exists an R with R > sRq(lo), R 2 > Tu 2 e~ 2 and u/e > 
Setting P = T l l 2 u/e, we need uT 1 / 2 > e 2 Ra(ui), and > e 1 " K T K / 2 , so 
these bounds hold for all sufficiently small e. □ 
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5. Invariance principle. In this section we prove the invariance principle 
Theorem 1.1. We assume that the conductances \i e are defined on the space 
(Q,P) where 

fi= [l,oo] Bd . 

We write /i e (w) = w(e) for the coordinate maps, and make the following 
assumptions on the environment (/U e ). 

Assumption 5.1. (1) (p e ) is stationary, ergodic, and invariant under 
symmetries of Z d . 

(2) n e G [1, oo) for all e G E d , P-a.s. 

(3) The conclusions of Theorem 1.2 hold for the VSRW associated with 

(Me)- 

As explained in the Introduction, our basic approach is to construct the 
"corrector" x '■ ^ x ^ d — * so that, for P-a.a. oo the discrete time process 

(5.1) M n = X n - X (u,X n ) 

is a P^-martingale with respect to the filtration T n = a(Xk,0 < k < n). 
The key steps in the proof of the invariance principle are: 

1. Tightness (a consequence of Theorem 3.5); 

2. The invariance principle for the martingale part. This is standard and 
follows from the ergodicity of our environment (see [29] , proof of Theorem 
2.6); 

3. The almost sure control of the corrector, for which we use the the er- 
godicity of the environment, the properties (4.11) and (4.12) and the 
quenched heat kernel estimates in Theorem 1.2 (see [34], or [12], Theo- 
rem 2.3). Note that all we need here is the ergodicity of the environment; 
ergodicity under the action of each direction as stated in [34], Remark 
1.3, is not required since one can use the cocycle property of the corrector 
(see [13, 26]). 

We now give the details. Let 

Qq = {uj: u(e) G [1, oo) for all e}. 

Since u>(e) satisfies Assumption 5.1 we have P(Oo) = 1. We write uj = (oj(e),e G 
Ed), and u(x, y) = oj({x, y}). For x G 7L d define T x : Q, — > by 

T x (u)(z, w) = lj(z + x, w + x). 

Let X be the VSRW with generator Ly given by (1.3), and qf(x,y) be the 
transition density of X . As v x = 1 is the invariant measure for X, 

<#(x,y) = PZ(X t = y) = tf(y,x). 
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Write 

(5.2) Q xy (u)=q^(x,y), QW( W ) = q%(x,y), x,yeZ d ; 
and note that Q xy < 1 for all x,y, with ^2 y Q X y = 1- We have 

(^•3) Qxy °T Z = Qas+z,y+zi Qxy = Qyx ■ 

We define the process Z, which gives the "environment seen from the parti- 
cle," by 

(5.4) Z t = T Xt uj, t€[0,oo), 

and define the discrete time process Z by Z n = Z n , n G Z+. 

Let W = L p (Q,F). For F € L 2 write F x = F oT x . Then Z has generator 

LF(u) = ^ Qox(u)(F x (lo) - 

Set 

5 (F, G) = E Q te (F - F X )(G - G x ). 

xez d 

Lemma 5.2. We have for F G L 1 , 

(5.5) EF = EF X ., 

(5.6) E(Q xi ? x)=E(Qo,-xi ? ). 

Proof. Since P is invariant by T x the first relation is immediate. As 
(Qo,*)-* = Q-x,o = Qo,- x by (5.3), E(Q 0x F x ) = E((Q 0x )_ x F) = E(Q ,- X F), 
proving (5.5). □ 

Lemma 5.3. If F,G € L 2 , then £(F,F) < oo, £(F,G) is defined, and 
LF G L 2 . 

Proof. Let FgL 2 . Then 

£(F,F)=E Qox{F-F x ) 2 

xez d 

<2EY,Qox(F 2 + F 2 ) 

x&Z d 

= 2EF 2 + 2E Y Qo* F x 

x&Z d 

= 2EF 2 + 2E ^ Qo- x F 2 = 4||F||§. 

x€Z d 
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Hence £(F, G) is denned for F, G G L 2 . Also, if F G L 2 , 

E\LF\ 2 =Ej2QoxQoy(F x - F)(F y - F) 

x,y 

1/2 / \ 1/2- 



< E 



^2 QoxQo y (F x — F) 2 j ( y^,QoxQoy(F y — F) 2 \ 



x,y ' x x,y 

2 



= £(F,F)<A\\F\\l □ 

Lemma 5.4. Let F,GeL 2 . Then 

(5.7) E(GLF) = —£(F,G). 

Proof. Using (5.5) we have 

(5.8) E(Qo,- x G(F_ x - F)) = E(Q 0x G x (F - F x )). 
So 

E(GLF) = Y,®GQ 0x (F x -F) 

x&L d 

x€Z d x£l d 

= \J2 E Qo*(GF x -GF + G X F - G X F X ) = —£(F, G), 

xez d 

where we used (5.8) in the last line. □ 



Now we look at "vector fields." We define for G = G(u, ij^xZ^K, 

EG = Y,^QoxG(-,x). 

x 

Definition. We say G(oj,x) has the cocycle property (see [13, 26]) if 

(5.9) G(T x u,y-x) = G(u),y) -G(u,x), P-a.s. 

Let H = L? be the set of vector fields G with the cocycle property and 
||G|| 2 =EG 2 <oo. 

Lemma 5.5. Let G = G(u,x) G L 2 . 

(a) G(lo, 0)=0 ; and G(T x uj,—x) = —G(uj,x). 

(b) If xo, x\, . . . , x n G 7L d , then 

(5.10) y~) G{T Xi _ x oj, Xj - Xi-i) = G(uj,x n ) - G(uj,x ). 
i=i 



INVARIANCE PRINCIPLE FOR THE RANDOM CONDUCTANCE MODEL 33 

Proof, (a) follows immediately from the definition. For (b), as G has 
the cocycle property 

G{T Xi _ x u, Xi - Xi-x) = G((jj, Xi) - G(uj, 

giving (5.10). □ 

It is easy to check the following: 

Lemma 5.6. 1? is a Hilbert space. 

For F £ L 2 we set 

VF(w, x) = F{T x uj) - F(u). 

Lemma 5.7. If F e L 2 , then VF eL 2 . 

PROOF. First, 

I|VF| 2 =J2^Qox(F x - F) 2 = 8(F, F) < oo. 

X 

Also, 

VF(T x u, y-x)= F(T y _ x T x uj) - F{T x u) 

= F(T v oj) - F{T x uj) = VF(w, y) - VF(u, x), 
so VF has the cocycle property. □ 

Lemma 5.8. Let G £L 2 . Then 

(5.11) eJ2Q^G(u,x) 2 <n\\G\\l 

X 

Proof. Write a 2 for the left-hand side of (5.11). Then using (5.9), 

(5.12) a 2 n = E ]T Q&^QxyiGincu, y - x) + G(u>, x)f . 

x y 

We now expand the final square in (5.12) and compute the three terms 
separately. We have 

E E E (u)Qxy(")G(T x GJ, y-x) 2 

x y 

= E E E Q-*,o ( T ^)Qo, y -x (T x uj)G(T x cj, y - x) 2 

x y 
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= E E E Q-*,o («)Qo»G(«, z? 

X z 

= eY,QoAu)G(oj,z) 2 = \\G\\ 2 . 

z 

Also, 

E E E ^ x) 2 = E £ Q£- 1} ( W )G( W> x) 2 = a 2 _! . 

x y x 

Finally, 

E E E (w)Q*»(o;)G(w, x)G(2> 

S J/ 

= E EE ^0x _1) (w)Qo,.*(2>)G(u;, x)G(T x u;, z) 

a; z 

< ('E^^Q(r 1) HQo, z (T,a;)G( W ,x) 2 ] 



1/2 



/ \ 1/2 

x fE^^QSr 1) HQo, 2 (T^)G(T^,z) 2 J 

^ X z ' 

= a n -i\\G\\. 

Thus a n < a n _i + ||G||, and so a n < ra||G||. □ 

Note that the following lemma uses the heat kernel lower bounds. 

Lemma 5.9. Let G € L 2 and 1 < p < 2. Then there exists a constant 
c p < oo smc/i i/iai 

(5.13) (E|G(-,x)n 1 /p<( Cp | x |)|| G ||. 
It follows that, F-a.s., 

(5.14) lim max |G( ^' x)l = 0. 

»woo| x |<„ n d+i 

Proof. By (5.9) and the triangle inequality we have 

(E|G(-,x)n 1 ^<|x|(E|G(-,e 1 )n 1 ^; 

so it is enough to bound E|G(-,ei)| p . By Theorem 1.2 there exists an in- 
teger valued random variable Wq with Wq > 1 such that P(Wq = n) < 
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c\ exp(— C2n s ) for some 5 > and qf(0,x) > cj,t~ d l 2 for t > Wo. Write £ n = 
g^(0,ei). Then 

oo 

(5.15) E\G(; ei )\P = J2nG(; ei )\n {Wo=n) . 

n=l 

Let a = 2/p, and let a' = 2/(2 — p) be its conjugate index. Then using 
Holder's inequality and (5.11), 

E\G(; ei )\n {Wo=n) 

= E(e n /a \G(;e 1 )n~ 1 / a l (Wo=n) ) 
<(E^G(-, ei ) 2 ) 1 /«(EC Q ' /Q l {iyo=n) ) 1 / a ' 

y ' 

< {n\\G\\ 2 ) 1/a c 4 n d/2a exp(-c 5 n 5 ) 

= c A n ( - d+2 ^ 2a exp(-c 5 n s )\\G\\ p . 

Summing the series in n we obtain (5.13). 
Using (5.13) with p = 1 we have 

F(max\G(u},x)\ > X n ) < (2n) d maxF(\G(u,x)\ > A n ) 

<cn d A" 1 maxE|G(a;,x)| < cn^A" 1 ||G||. 

|a;|<n 

Taking A n = n rf+3 and using Borel-Cantelli gives (5.14). □ 

Following [29] we introduce an orthogonal decomposition of the space L 2 . 
Set 

L 2 p = cl{VF,FeL 2 } inn, 

and let L 2 be the orthogonal complement of L 2 in %. (Here p stands for 
"potential" and s for "solenoidal.") 

Lemma 5.10. Let G G L 2 . Then for each x, EG(x,w) = 0. 

Proof. Fix x G Note first that if G = VF, where F G L 2 , then 
EG(w, s) = E(F C _- F) = EF-c - EF = 0. 

Now let G £ L 2 . Then there exist F n G L 2 such that G = lim n VF n in 
L 2 . Since P(Qox > 0) = 1, it follows that VF n (uj,x) converges to G(u,x) 
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in P-probability. By Lemma 5.9, for each p£ [1,2) the sequence X7F n (uj,x) 
is bounded in L p (f},P), and therefore VF n {u,x) converges to G(u,x) in 
L 1 (S7,P). So EG(u,x) = lim n EVF n (uj,x) = 0. □ 

We define the semi-direct product measure P* = P x P® . 

Lemma 5.11. LetGeL%. Then 

(5.16) ^ Q 0x (uj)G(u,x) = 0, P-a.s. 

x&Z d 

Hence M n = G(oj,X n ) is a P® -martingale for P-a.a. u. Further, writing 
\\G(u;,-)\\ 2 = ^Q 0)X (uj,x)G(u;,x) 2 , 

we have 

n— 1 n—1 

(5.17) (M) n = ]T \\G(T^u, -)H 2 = E H G &. Of- 

k=0 k=0 

Hence 

(5.18) E*(M n ) 2 = n||G|| 2 . 

Proof. If F e L 2 and G e I 2 , then using Lemma 5.5, 

£ EQ 0a; G(w,x)F x = £ EQ 0;E (T^)G(T_^,x)F x .(T._ x .w) 



EEQo,-xH(-G(a;,-x))Fl 
- ^2 EQ Qx (uj)G{uj,x)F(uj). 



x& 



Thus 

(5.19) J2 E QoxG(-,x)(F + F x )=0. 

x&Z d 

If Gel 2 , then 

= E(GVF) = E ^QoxG(; x){F x - F), 

x 

and so K^2Qq x GF = 0. Since this holds for any F £ L 2 we obtain (5.16). 
To show that M is a martingale it is enough to prove that for any x, 

(5.20) E°(G(Lu,X n+1 )-G(LO,X n )\X n = x) = 0. 
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However, using (5.16), 

E%(G(u,X n+1 )-G(u,X n )\X n =x) 

= ^2Qxy{u)(G(LU,y) - G(u),x)) 
y 

= ^ Qo,y-x(T x u;)G(T x u), y-x)=0. 
y 

Recall that (M) is the unique predictable process so that M 2 — {M) n is a 
martingale. We have 

E'*(M 2 +1 - M 2 \X n = y) = EZ((M n+1 - M n ) 2 \X n = y) 

= Y,Qyz(u)(G(u,z)-G(u,y)) 2 

z 

= ^2Qo,z~y(T y uj)(G(z -y,u)) 2 

z 

= \\G(T y co,-)\\ 2 , 

and (5.17) follows. 
Finally, 

n-l 

E*Af n 2 = E(i?X) = E(E°jM) n ) = £ n\G(T2 k u>, -)f = n\\G\\ 2 . □ 

fc=0 

Let II : Mr — > M d be the identity, and write LL,- for the jth coordinate of IT. 
Then Hj(y — x) = ILj(y) — Hj(x), so ILj has the cocycle property. Further by 
(4.14), 

I|IL;| 2 = E^2Q 0x \xj\ 2 <oo, 

X 

so IL; £H. So we can define Xj G L 2 and 3> j £ L 2 by 

V-j = Xj + $j£Ll®L 2 s ; 

this gives our definition of the corrector \ = (xi> • • ■ i Xd) : ^ x Z d — >■ M d . We 
will sometimes write x(;c) f° r xC'j^O- Note that conventions about the sign 
of the corrector differ — compare [34] and [12]. As the environment process 
is invariant under isometries of Z d , ||$j|| = ||$i|| for each j = 1, . . . ,d. 
The following proposition summarizes the properties of x an d 3>- 

Proposition 5.12. (a) M n = X n - x(u, X n ) is a -martingale. 
(b) For each x £ Z d , x(-, x ) G L 1 . 
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(c) For each j = 1, . . . , d 

E^Q xH|^(w,x)| 2 = ||$i|| 2 < oo. 

x 

(d) x i- s sublinear on average; for each e > 

(5.21) lim n~ d V h\ x (u tX )\ >en ) = 0, F-a.s. 

\x\<n 

Proof, (a) and (b) are immediate from Lemmas 5.11 and 5.9, and (c) 
is immediate from the definition of 3>j as a projection in 1? . Let e\ be the 
unit vector e\ = (1, 0, . . . , 0). By Lemma 5.10 we have Ex(-, ei) = 0. So since 

n 

(5-22) x{u,nex) = ^x(T (A ._ 1)ei o;, ei) 

fc=i 

and as x has the cocycle property, the ergodic theorem implies that lim n n _1 x(w, 
nei) = P-a.s., and (d) then follows by the results in Section 6 of [26]. □ 

Lemma 5.13. The processes Z and Z are ergodic under the time shift 
on the environment space fi. 

Proof. This is well known; see [17], Lemma 4.9, and Section 3 of [10] 
for a careful proof in discrete time. □ 

Proof of Theorem 1.1. We begin with the VSRW. The arguments 
are very similar to those in [10, 12, 34], so we only mention the key points. 
We define 

(5.23) M n = $(u,X n ), M t (e) =£M Li/£ 2j, t>0, 
so that 

(5.24) xl e) =eX Vt/e ^ =M} e) +e X (u;,e- 1 xi E) ). 

Thus it is sufficient to prove that the martingale converges to a mul- 
tiple of Brownian motion, and that for P-a.a. cj, the second term in (5.24) 
converges in P^-probability to zero. 

We start with the control of the corrector, and use [12], Theorem 2.4. 
This proves that if the corrector \ has polynomial growth, and is sublinear 
on average, then Gaussian upper bounds on the heat kernel imply pointwise 
sublinearity of x- Thus, using (1.10), (5.14) and (5.21) we have that for 
P-a.a. u, 

(5.25) lim max lx( ^ ,x)l = 0. 

n-s>oo| x |< n n 
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Given (5.25) the Gaussian upper bounds then imply that, for P-a.a. lj, 

(5.26) £x(oj,Xu/ £ 2\) — >■ in P°-probability. 

For the convergence of M( £ \ we proceed as in [10]. Let v G R d be a unit 
vector, write M% = v ■ M n , and let 

F K (u) = I%(\Mt\ 2 -,\M?\>K). 

Then Fk is decreasing in K and 

EF K <EF < dW^xW 2 . 

In the notation of Lemma 5.11, Fq(u) = \\v ■ $>(oj, -)\\ 2 , and so by (5.17) the 
covariance process of M v is 

n-l 

(M*) n = ^F (Z k ). 

k=0 

So by Lemma 5.13 we have n~ 1 (M v ) n —^MFq, P^-a.s., for P-a.a. u. 

Using the same arguments as in [10], Theorem 6.2, it is straightforward 
to check the conditions of the Lindeberg-Feller FCLT for martingales (see, 
e.g., [20], Theorem 7.7.3), and deduce that v ■ converges to a constant 
multiple of Brownian motion. Hence M^ £ ' converges to an Revalued Brow- 
nian motion with nonrandom covariance matrix D given by Dij = E<I>j < l > j. 
Since the law of the random variables oj(e) is invariant under symmetries of 
Z d , we deduce that there exists <jy > such that D = <JyI, and that 

(5.27) cr^=I<^. 

This establishes the convergence of X^; using Lemma 4.12 gives the con- 
vergence of X^ to the same limit. 

The global upper bounds on qf(0,x) in Lemma 2.8 imply that if A > 
and \t 1 / 2 > 1, then 

P£(\Xt\ < At 1/2 ) < ct- d / 2 \B{0, At 1 / 2 )] < c\ d . 

Hence there exists A > such that for all large t, 

P»{\X t \>\t l ' 2 )>\, 

which implies that <iy > 0. 

We now consider the CSRW. Recall from (1.1) the definition of /J, x (u), set 
F(ui) = Ho(uj), and 

(5.28) A t = I ^x s ds= I F{Z s )ds. 

Jo Jo 
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Then if n = inf{s > : A s > t} is the inverse of A, the time changed process 

(5.29) Y t = X Tt 

is the CSRW. 

By the ergodic theorem for the process Z, 

lim t~ x A t = EF = 2 cffiu e , P*-a.s. 

t— >oo 

So if Eu e < oo then Tt/t—> a a.s. where a = 1/2 cffiu e > 0. Let Y^ = eY t j e i . 
Then 

(5.30) Y^=xi? + (Y t V-X$) 

and using Theorem 4.11 we have for any fixed to > that 

(5.31) sup \Y t {£) -X { a l>\ 

o<t<t 

converges in P® -probability to 0, for P-a.a. to. Thus Y^ converges to (TcB' t 
where B' is a Brownian motion and = aoy > 0. 

In the case when E/u e = oo we have that Tt/t — > 0, and hence Y& converges 
to a degenerate limit. □ 

We conclude this section by stating a local limit theorem for qt(x,y) (for 
the VSRW). Write 

hix) = (27rta 2 v )- d / 2 e-^ 2 ^v t 

for the Gaussian heat kernel with diffusion constant <jy where <jy is as in 
Theorem 1.1. 

Theorem 5.14. Let X be the VSRW. Let T > 0. For x £ R d write [x\ = 
([xi\,. ■ ■ , [x d \). Then 

(5.32) lim sup sup|n d/2 ^(0, [n l/2 x\) - k t (x)\ = 0, F-a.s. 

Proof. This is proved as in Section 4 of [3]. We have to verify As- 
sumptions 4.1 and 4.4 in [3], but this is straightforward given the invariance 
principle and heat kernel bounds in Theorems 1.1 and 1.2, and the PHI in 
Theorem 4.7. 

Note that as v is the invariant measure for X, in Assumption 4.1(d) all 
we need is that v(A n (x,r)) j \2n l l 2 r)d converges, and as v x = 1 for all x; this 
is easy. (Here A n (x, r) = (in 1 / 2 + [-rn 1 / 2 , rn 1 / 2 ] d ) n Z d .) □ 
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Remark 5.15. In this section we have constructed a corrector x(w,x) 
so that the process 

(5.33) M n = X n - X {u,X n ), neZ + , 

is a (discrete time) martingale. It is natural to ask if x( w i ") is a ls° a corrector 
for the continuous time process Xf. 

For the RCM with i.i.d. conductances it is straightforward to check that 
the condition in Lemma 2.21 involving the set A{K) holds P-a.s. (see [12], 
Lemma 3.1, for a similar argument). We can therefore use Lemma 2.21 with 
M') = xi 00 , ") t° deduce that 

(5.34) M t = X t - X (tJ,X t ), t £ M + , 
is, for P-a.a. oj, a P^-martingale. 

6. General ergodic environments. We conclude this paper with some 
remarks on more general ergodic random environments. First, note that the 
proof of the invariance principle in Section 5 just uses the facts that the 
environment is stationary, symmetric and ergodic, and that the heat kernel 
bounds in Theorem 1.2 hold. 

In the proof of Theorem 1.2 the full strength of the assumption that [i e 
were i.i.d. was only used at one point, in Theorem 4.3, where we controlled 
the probability that a ball was not very good. The heat kernel upper bounds 
in Section 2 only require Assumptions 2.1 and 2.6, together with a compar- 
ison of the metrics d(x,y) and d(x,y). Given these upper bounds, and using 
the fact that /x e > 1, no additional hypotheses on fj, e were needed to obtain 
the lower bounds in Section 3. We therefore have the following: 

Theorem 6.1. Let fj, e , eGEj be a stationary symmetric ergodic envi- 
ronment, satisfying for some c\ > 0, 

(6.1) jU e £[ci,oo) for all e& E^F-a.s. 

Let d u {x,y) be the metric given by the first passage percolation construction 
of (2.6), and (as in Definition 2.9) let V X (X) be the smallest integer such that 
(x,V x (X)) is X-very good. Suppose that there exists Xq < oo and rj £ (0,1) 
such that 

(6.2) nV x (X Q )>n)< Cl e-^ n \ 

Then the conclusions of Theorems 1.1, 1.5(a)— (c), 1.3, 4-7 and 5.14 all hold 
for the environment (/x e ). 

Proof. We begin by considering the heat kernel bounds in Theorem 
1.2. As in Lemma 4.1, it is immediate that Assumptions 2.1 and 2.6 hold for 
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(jU e )j F-a.s. Using the hypothesis (6.2) instead of Theorem 4.3, the arguments 
in Section 4 [except for Theorem 4.9(a), for which see Remark 6.2 below] 
hold in this more general context, and give Theorems 1.2 and 4.7. 

Given Theorem 1.2, the arguments in Section 5 then give the invariance 
principle (Theorem 1.1) and local limit theorem (Theorem 5.14). 

Combining these results gives the Green function estimates in Theorem 
1.3. □ 

Remark 6.2. The proof of Theorem 4.9 used the fact that the bounds 
in Theorem 4.6 hold for 1 — 77 as well as for 77. If we only have (6.2) then we 
obtain 

(6.3) E<#{x,y)<c 1 t- d / 2 e- C2 \ x ~ y \ 2 / t if t > c 3 V \x - y\ 1+r > , 

(6.4) ^{x,y)>c A r dl2 e- c ^ x - y \ 2/t if t > c 6 V \x - y\ 2 ~ r > . 

Remark 6.3. If fx e is bounded and bounded away from 0, so there exist 
< c\ < C2 < 00 such that P(/i e £ [ci,C2]) = 1, then the metrics d(x,y) and 
d(x,y) are comparable. So, taking Ao large enough, (6.2) holds. 

Remark 6.4. See [12], Lemma 3.1, or [30], Lemma 5.3, for percolation 
arguments which are more robust than Theorem 4.3 and which may be useful 
for establishing (6.2) in more general contexts. 

Remark 6.5. If fj, e is stationary and ergodic, but not invariant with 
respect to symmetries of Z d , then if (6.2) holds, we still obtain Theorem 1.2, 
and the convergence of to a Brownian motion with covariance matrix 
D. However, D need not be diagonal. 

Remark 6.6. Unlike ergodic bounded conductance models, the results 
of this paper certainly do not hold for all unbounded stationary symmetric 
ergodic random environments. For example, let cf = 2,3,4 and let T be a 
uniform spanning tree on Z rf (see [9]). Then T is 1-sided, so from each x £ 7j d 
there is a unique self-avoiding path 7 X to infinity. Let a(x) be the first point 
on this path. Then a : 7L d — > 7L d and the path 7^ is {x, a(x), a 2 (x), . . .}. 

Let N(x) be the set of points in T which are disconnected from infinity 
by deleting the bond {x,a(x)}, and let n(x) = |iV(x)|. As x £ N(x), n(x) > 1 
for all x. Let \x e = 1 for edges e £ which are not in T ■ Each edge e £ T is 
of the form e = {x,a(x)} for some x, set 

V{x,a(x)} = n{x)e n{x)2 . 
Let Tj, % > 1, be the jump times of the VSRW X . Then 



(6.5) P*(X Tl ^a{x)) 
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Fix x, and let the neighbors of x in T be a(x),yi, ■ ■■ ,yk- Then 

k 

^ y = {2d-k-l) + Y,n{yi)e n ^ 2 . 

yj^a(x) i=l 

Since fi Xja (x) = n(x)e n ^ 2 and n(x) = 1 + X) n (2/i)> it is easy to see that 
Po ( x ) = P*(X Tl ^ a(x)) < 2de~ nix)2 +maxe' ife)2 - n(l)2 < 2de~ n{x)2 + e~ n{x) 1 d . 

i 

So ^kPoi^ix)) < oo, and it follows that ultimately the process X moves to 
infinity along a path j x for some x. Since Y^kf i a k (x) a k + 1 (x) ^ 00 ^°* s takes 
finite time. Hence the quenched invariance principle Theorem 1.1 fails, as 
well as the Gaussian bounds in Theorem 1.2. 
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